Statistical theory of hierarchical avalanche ensemble 



Alexander I. Olemskoi 
Physical Electronics Department, Sumy State University 
2, Rimskii-Korsakov St., 244007 Sumy UKRAINE 
q^" E-mails: Olemskoi@ssu.sumy.ua; Alexander@olem.sumy.ua 

ON ' The statistical ensemble of avalanche intensities is considered to investigate diffusion in ultra- 

metric space of hierarchically subordinated avalanches. The stationary intensity distribution and 
the steady-state current are obtained. The critical avalanche intensity needed to initiate the global 
avalanche formation is calculated depending on noise intensity. The large time asymptotic for the 
probability of the global avalanche appearance is derived. 

Q . I. INTRODUCTION 

In recent years considerable study has been given to the theory of self-organized criticality (SOC) that explains 
avalanche dynamics in vast variety of systems type of sandpile model 0, intermittency in biological evolution ||, 
earthquakes and propagation of forest-fires, depinning transitions in random medium and so on (see Without 
any special tuning of their external parameters, such systems evolve to a critical state that displays universality in 
both spatiotemporal behaviour and avalanche distribution over size s and intensity / . So, the latter is given by 
homogeneous function 



C3 



T3 

O 
o 



> 



c3 



P(s,f) = s- T g(x), x = s/s c , (1) 

where a cutoff size s c is determined as 

s c ~c/wr 1/ff , (2) 

f c is a critical intensity, r, a are critical exponents. Such relations are obtained as result of making use of scaling- 
type arguments supplemented with extensive computer simulations 0] as well as within framework of the theory of 
branching processes H. These methods allow to find the magnitudes of the parameters f c , t, a, but as to the function 



g(x) in Eq.(|l]), it is known only that this function is a monotonically falling down 
, This work is devoted to determination of the scaled function g(x) within the framework of statistical theory that 
deals with avalanche ensemble in the course of SOC progressing. The cornerstone of our approach is hierarchical rep- 
ON \ resentation of the avalanche ensemble that is embedded into ultrametric space corresponding to hierarchical landscape 
of the system energy Q. 

The paper is organized as follows. Sec.|| deals with the theory of hierarchical coupling between elementary 
avalanches arising from the complexity of the phase space landscape of self-organized system. The time dependence 
of the probability of the global avalanche appearance is studied in Sec. [II. Sec.|^ contains discussion of obtained 
results. 
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II. INTENSITY DISTRIBUTION IN HIERARCHICAL AVALANCHE ENSEMBLE 

Let the maximum number of hierarchically subordinated avalanches N be on the bottom hierarchical level s = 
rS | where the avalanche intensity equals /. Correspondingly, there is the only avalanche of the intensity F > / on the 
top level s — s c ^> 1. The problem is to find the dependencies N(s) and F(s) that define the distribution of avalanche 
number and their intensity over hierarchical levels s £ [0, s c — 1]. 

The first part of the problem can be approached in terms of geometry by representing the avalanche ensemble as 
a hierarchical Cayley tree fi],|||. The basic types of the trees are shown in Fig.l: regular tree with integer branching 
ratio j, regular Fibbonachi tree with fractional one j = r rj 1.618, degenerate tree with the only branching node per 
level and the tree of our primary concern - irregular tree. Let k be the numbering index for the levels, so that k 
increases from the top level to the bottom one. The variable 

s = s c — k (3) 
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then defines the distance in the ultrametric space ||. Geometrically, objects of this space correspond to the nodes of 
the bottom level (k — s c ) of a Cayley tree. Since the distance between the nodes is defined by the number of steps to 
a common ancestor, the distance is eventually the level number (||), counted from the bottom. 

As it can be seen in Fig. la, in the simplest case of regular tree with integer branching ratio j the number of 
avalanches Nk — j exponentially decays to zero with the distance s between them: 

N(s) =Nexp(-slnj), N=f c . (4) 

In Eq.(^) the equality (|^) is used and the avalanche number N is related to the total number of levels s c . For the 
Fibbonachi tree (see Fig. lb), where Nk = qr k , q w 1.171, r w 1.618 ||, we have 

N(s) = Nexp(-s\nr), N = qT s <=. (5) 

When Eq. (||) is compared with Eq. ([|), it is clear that the exponential decay remains unaltered in the case of 
fractional branching ratio and characterizes the regularity of tree. 

For the degenerate tree (see Fig.lc) Nk = (j — l)k + 1 and Eq.(|^) provides the following linear dependence 

N(s) = N-(j- 1) 8 , N=(j- l)s c + 1. (6) 

It can be shown that in the case of irregular tree, displayed in Fig. Id, the power law dependence is realized: 

N k = fc a , a > 1. (7) 

Indeed, the latter can be regarded as an intermediate case between the exponential Eqs.(^), (|J) and linear Eq.(|^) 
obtained for the limiting cases of regular and degenerate trees, respectively. Formally, the approximation means 
that a function N(x) defined on the self-similar set of hierarchically subordinated avalanches is homogeneous, N(kx) = 
k a N(x). It is convenient to rewrite Eq.(^) in term of the distance: 

N k = N(l - s/s c ) a , N = s a c , a>l. (8) 

Now, let us define Fk as total intensity of avalanches on the fc-th level. If one considers the value Fk as an effective 
density of some particles and the level number A; as a coordinate, then corresponding density of hierarchical current 
can be taken in the Onsager-type form: 

3k = -D(Fk) (9) 

Here, within the multiplicative noise approach, the effective diffusion coefficient has form fic| l 

D(F) = DF~ a (10) 

to depend on the constant D > and the exponent a. 

The basic assumption of this section is that the total current J of all avalanches at given hierarchical level is 
independent on the level number k: 

jkNk = const = J. (11) 
Inserting Eqs.(^|)-(^0|) into Eq.([ll|) gives the total avalanches intensity on the level k 

F k = Fk-\ ^-=1^4 b^^-l>0 (12) 
a — 1 L> 1 — a 

normalized by the maximum value F = Fk=\. Introducing the distance (||), we obtain 

F(s) = f(l-s/ Sc )-\ (13) 
where the intensity at the bottom level s — is 

/ = Fs- b = FN- b/a . (14) 
After generalizing Eqs.(ft2h, (|14|), the following scaling relation can be assumed 
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F k = N b ' a k- b f k (15) 
where fu is a slowly varying function. According to Eqs.(|9|)-(|ll|), this function obeys the Landau-Khalatnikov equation: 



where one denotes 



= \nk b 



and the effective potential is introduced 
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fl~ a = (J/bD)N-( a -W" 
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(16) 
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As indicated in Fig. 2, the potential V reaches its maximum value V c = (1 — a)/2(l + a) at x = 1 and decreases 
indefinitely at £ > 1. So, in order to initiate the global avalanche formation, a low intensity avalanche with f < f c 
at the bottom level needs to penetrate the barrier V c of the potential (|l8|). It implies fluctuation mechanism for the 
SOC regime progressing provided that a; is a stochastic variable and we proceed with Langevin-type equation derived 
from Eq.(|l6|) by adding a Gaussian white noise to the right-hand side: 



da; _ dV + 
dft dx ' 



(19) 



(0=0, (CWC(«')) = 22W(k - *'). 



(20) 



where the noise intensity D equals the diffusion coefficient in Eq. (|i0|) . 

The usual way to study a set of solutions to the stochastic equation ( |l9|) is to introduce distribution function g(n, x) 
associated with the probability of solution's realization. It is known that g(n,x) obeys the Fokker-Planck equation 
& 



dn dx 



0. 
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dV dg 
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(21) 



Since there is no current at the equilibrium state (j 
intensities at the bottom level 



0), the corresponding distribution function of avalanche 



g (x) oc exp(-V(x)/D) 



(22) 



is dictated by the potential (|18|) . In the case of non-equilibrium steady state the probability density g does not depend 
on the hierarchical level variable k and the current j being constant, in compliance with the conservation law jll]), 
can take a non-zeroth value. From Eq.(21) the stationary distribution then is expressed in terms of the equilibrium 
distribution go(x) and the current j fTT 
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where the boundary condition g — > as / — > oo is taken into account. 

Given the intensity / equation ( p3| ) allows the current j to be found. In trying to do it, special consideration should 
be given to the fact that the intensity / is bounded from below, / > G, by a gap G that is inherent in hierarchical 
ensemble of avalanches [[|. Indeed, after merging of avalanches within a hierarchical cluster of the size s g , all s, 
such that s < s g , are appeared to be dropped out the consideration as well as low intensities with / < f(s g ) = G 

with the second boundary condition 



(see Fig.l). The expression for the current j then can be derived from Eq.(|2 
g(G) = go(G). The result reads 



j = 2DW { 1 + erf 



1-a 



2D V f c 



(24) 
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where the factor 



1 1 - a 

W cx expt-KAD), V c = - — — (25) 

z 1 + a 



gives the probability that fluctuation will surmount the barrier V c of the potential (|18|). Equation (|24|) shows that 
in the case of small gap, G <C f c , the current j is about WD, but the current is doubled under G — f c . It can 
be understood if we picture the effect of the gap as a mirror that reflects diffusing particles at the point / = G: if 
G « J c a particle penetrating the barrier can move along both directions, but in the case of G = f c the mirror is 
placed at the point corresponding to the top of the barrier and all particles go down the side where the intensity / 
grows indefinitely. 

Given the current j the stationary distribution function <?(/) is defined by Eq.([2^), according to which, g(f) ~ <?o(/) 
in the subcritical region / < / c , while in the supercritical range / 3> f c we have go{f) 3> <?(/) due to indefinite increase 
of go(f). As far as the stationary distribution is concerned, it can be derived from the current definition (pl|), where 
the last diffusion term is negligible for supercritical intensities: j w — (dV/dx)g. The result is that the probability 
g(f) remains almost unaltered (g(f) ~ $(/c)) in the range from the critical value f c up to the boundary one f g and 
#(/) at / 3> f g (sec [^2|). The growth of f g is governed by the equation 

it =Dk jr (26) 

that is a counterpart of the known gap-equation presenting the system behaviour below the critical value f c . 

Since the above picture is essentially statistical, it enables the critical avalanche intensity f c for the transition point 
to be found. Indeed, when the definition of the macroscopic current J in Eq.(pd|) is compared to that of the microscopic 
current j in Eq.(^lj), it is apparent that they differ from one another only by the factor 7V( a ~ 1 '/ a = dependent 
on the total number of avalanches N. On this basis, the last expression of Eq.(|l7|) and Eq.(p4|) at G — 0, D <C 1 give 
the desired result: 

f c = Fexp(-D /D), D« 1 = 2{\ + a), (27) 

where the pre-exponent factor F determines the probability of the barrier penetrating and cannot be calculated 
within the framework of the presented approach. Equation (|27j) bears a resemblance to the well-known result of the 
superconductivity BCS theory for the temperature of the phase transition and predicts the slow growth of the critical 
intensity f c of elementary avalanche with the hierarchical diffusion coefficient D that plays the role of the parameter 
of effective interaction. 



III. KINETICS OF THE GLOBAL AVALANCHE FORMATION 



Since the ensemble of hierarchically subordinated avalanches represents a self-similar set, the probability distribution 
P(s, /) of avalanches in the course of SOC process is a homogeneous function ([j]), where g(s/s c (f)) = g(f) is the 
stationary distribution of elementary avalanches considered in the previous section. Physically, Eq. (0) implies that the 
total intensity F, being measured by the scale N h / a , equals the intensity of an elementary avalanche / in accordance 
with Eq.((l|). 

In this section we are aimed to describe kinetics of the global avalanche formation produced by virtue of the 
hierarchical coupling between elementary avalanches. As it has been clarified in Secjn], this process can be conceived 
of as diffusion in ultrametric space that makes the distribution (|l|) mounted. In order to find the conditional probability 
V{t) that no global avalanche will appear at time t one has to integrate over s the distribution (Q) weighted with the 
function 



Pa (t) = cxp (-</*(«)), *(«) = h exp (F(s)/D) 



(28) 



descriptive of Debay relaxation with the time t(s) governed by the barrier height F(s) (to is a microscopic time). By 
using the steepest descent method, it is not difficult to derive the late time (t —> oo) asymptotic formula 
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This equation has been obtained by assuming that the condition 1 <C s T , 
of the maximum of integrand and obeys the equation 



< s c is met, where s m denotes the location 
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Dr(l-y) 1 + b t ( f 



s 



exp l-j-(l-y)-' , y^— . (30) 
of y t \ D J s c 

Taking into consideration the scaling relation (|[) for the number of hierarchical levels s c , which is the cut-off parameter, 
we readily come to the conclusion that the condition is satisfied provided 

/-/c</, <»t e/ exp((/ c / J D)- 1/b -l) \ (31) 

So, the intensity / in Eqs.(]29|),(j30|) can be replaced by the critical value f c . In accordance with Eq. (p9|) , the probability 
V(t) = 1 — V(t) of the global avalanche appearance logarithmically increases in time up to the value V = 1 — (f c /D) T / b . 



In order for the probability V to be non-negative, the factor F in Eq.(27) must be equal Fq = (e/2)(l + a) 1 , whereas 



the effective diffusion coefficient D must be bounded from above by the value Dq = (1/2)(1 + a) 



IV. DISCUSSION 



According to the presented picture, the initiated elementary avalanches form statistical ensemble of hierarchically 
subordinated objects, characterized by intensity / and distance in ultrametric space s (the latter corresponds avalanche 
size [Q). Since the global avalanche formation is caused by effective diffusion in the space, then, similar to Brownian 
particle with coordinate / at time s, the ensemble can be described by Langevin equation (|l9| ) subjected to the noise 
Eq.(p0[) with D being the effective diffusion coefficient and corresponding to the Fokker-Planck equation (21). The 
stationary intensity distribution and the steady-state current are given by Eqs.(|23|), J24|). The condition of current 
conservation Eq.(ll) yields the avalanche intensity distribution ( |l3| ) over hierarchical clusters in the ultrametric space. 
The ensemble of elementary avalanches, being weakly dependent on s, is governed by the effective potential ( |l8|) that 
reaches its maximum at the critical intensity (^7|) (see Fig. 2). So, the global avalanche generation requires supercritical 
elementary avalanche intensity, / > / c , to surmount the barrier V c with the characteristic time (cf. Eq.(pq)) 



T «i o ex P 0/ c /£>), 



n 



Vc 2 1 + n 



(32) 



This picture bears some resemblance with the formation process of supercritical embryo in the theory of the first-order 
phase transitions (l2| , where in the course of transformation the next stage is the diffusion growth of the embryo. 
Analogously, in the case under consideration the above growth implies an increase of the supercritical avalanche in 
intensity F(s), Eq.(|l3|), due to the diffusion growth of hierarchical cluster in ultrametric space. As a result of the total 
cluster formation, we have the logarithmically slow large time asymptotic for the probability of the global avalanche 
appearance: 



V(t) = 1- V 
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t ef = {T/b)V 1/T t , 



(33) 



where time t is counted from the instant T, Eq.(|32|), and V is the maximum probability that no global avalanche will 
occur 
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(34) 



From Eq.(p4|) the probability is determined by the ratio of the noise intensity D (see Eq.(|2(j)) and its maximum value 
Do = (1/2)(1 + a) -1 . The key point is that the maximum probability V = 1 — V of the global avalanche appearance 
is completely suppressed under high intensity variance in ensemble of elementary avalanches (see Fig. 3). 
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FIGURE CAPTIONS 



FIG. 1. Different types of hierarchical trees (the level number is indicated at left, corresponding number of nodes - at 
right): a) regular tree with j — 2; b) Fibonacci tree; c) degenerate tree with j — 3; d) irregular tree. 



FIG. 2. The effective potential (18) as a function of f/f c at a = 0.1. 



FIG. 3. The dependence of the maximum probability V = 1 — V of the global avalanche appearance on the intensity variance 
D in ensemble of elementary avalanches. 
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